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GLOBAL STABILITY OF STEADY STATES IN THE CLASSICAL STEFAN PROBLEM 


MAHIR HADZIC AND STEVE SHKOLLER 


ABSTRACT. The classical one-phase Stefan problem (without surface tension) allows for a continuum of steady state solu¬ 
tions, given by an arbitrary (but sufficiently smooth) domain together with zero temperature. We prove global-in-time stability 
of such steady states, assuming a sufficient degree of smoothness on the initial domain, but without any a priori restriction 
on the convexity properties of the initial shape. This is an extension of our previous result da in which we studied nearly 
spherical shapes. 


1. Introduction 

1.1. The problem formulation. We consider the problem of global existence and asymptotic stability of classical 
solutions to the classical Stefan problem, which models the evolution of the time-dependent phase boundary between 
liquid and solid phases. The temperature p(t,x) of the liquid and the a priori unknown moving phase boundary T(f) 


must satisfy the following system of equations: 

Pt — Ap = 0 in 0(f); (1.1a) 

V(T(t)) = -d n p on T(f); (1.1b) 

p= 0 on F(f); (1.1c) 

p(0,-)=Po,fi(0) = n. (l.id) 


For each instant of time f € [0,T], O(f) is a time-dependent open subset of with d> 2, and 17(f) =<90(f) denotes 
the moving, time-dependent free-boundary. 

The heat equation (1 1.1 at models thermal diffusion in the bulk O(f) with thermal diffusivity set to 1. The boundary 
transport equation (ll.lbb states that each point on the moving boundary is transported with normal velocity equal to 
—d n p=—S7p-n, the normal derivative of p on 17(f). Here, n(f,■) denotes the outward pointing unit normal to l’(f), 
and V(T(f)) denotes the speed or the normal velocity of the hypersurface l’(f). The homogeneous Dirichlet boundary 
condition (ll.lcb is termed the classical Stefan condition and problem (11.1b is called the classical Stefan problem. It 
implies that the freezing of the liquid occurs at a constant temperature p = 0. Finally, in (ll.ldb we specify the initial 
temperature distribution po : Q —> R, as well as the initial geometry SI. Because the liquid phase fl(t) is characterized 
by the set {igR 1 * : p(x,t) > 0}, we shall consider initial data po > 0 in SI. Thanks to (11. lab , the parabolic Hopf 
lemma implies that d n p(t ) <0 on T(f) for f >0, so we impose the non-degeneracy condition (also known as the 
Rayleigh-Taylor sign condition in fluid mechanics [ 43] [45 i:47 1 i f4l[17]IT6l 1 : 

— d n po > A > 0 on T(0) (1.2) 

on our initial temperature distribution. Under the above assumptions, we proved in lf27l that (11.1b is locally well-posed. 

Steady states (u,T) of (11.1b consist of arbitrary domains with I £ (f and with temperature u = 0. The main goal 
of this paper is to prove global-in-time stability of such steady states, independent of any convexity assumptions. Our 
analysis employs high-order energy spaces, which are weighted by the normal derivative of the temperature along 
the moving boundary; we create a hybridized energy method , combining integrated quantities with pointwise methods 
via the Pucci extremal operators, which allow us to track the time-decay properties of the normal derivative of the 
temperature. This hybrid approach appears to be new, and is a natural extension of our previous work (28], which 
necessitated perturbations of spherical initial domains. 
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1.2. Notation. For any s > 0 and given functions f :Q—> R, <p: T —>M we set 

ll/IU = ll/ll^(0) and |^| s ^ f ||^||^ (r) . 


def . ■ def 

If i = then f,i = d x if is the partial derivative of / with respect to x l . Similarly, f,ij = 8 x id x jf, etc. For time- 

differentiation, ft = dtf. Furthermore, for a function we shall often write /(f) for /(t, -), and /(0) to mean 

f{0,x). The space of continuous functions on 12 is denoted by C'°(f2). For any given multi-index a = (ai ,... ,a d ) we 
set 

d s = d? 1 ...d° d . 


We also define the tangential gradient 8 by 8f = V/ — 8n fN, where N stands for the outward-pointing unit normal 
onto <917 and djyf = N ■ V/ is the normal derivative of /. By extending N smoothly into a neighborhood of T inside the 
interior of 17 we can define 8 on that neighborhood in the same way. We employ the following notational convention: 


8f = (8 1 f,...,8 d f), 8 s f^ (, 8? /,..., 8 a d d f), 

where a= (ai.... ,a d ) denotes a multi-index. The identity map on 17 is denoted by e(x) — x, while the identity matrix 
is denoted by Id. We use C to denote a universal (or generic) constant that may change from inequality to inequality. 
We write X < Y to denote X < CY. We use the notation P(s) to denote a generic non-zero real polynomial function 
of s 1 / 2 with non-negative coefficients of order at least 3: 

m 

P(s)=^2 c iS~5- , Ci> 0, raeN 0 . (1.3) 

i =o 

The Einstein summation convention is employed, indicating summation over repeated indices. 


1.3. The initial domain Q and the harmonic gauge. For our initial domain O we choose a simply connected domain 
ft. C R d , where the boundary Off will be denoted by F. We further assume, without loss of generality, that the origin 
is contained in ft, i.e. 0 (E ll. We transform the Stefan problem (11.11) set on the moving domain f1(f), to an equivalent 
problem on the fixed domain 17; to do so, we use a system of harmonic coordinates, also known as the harmonic gauge 
or Arbitrary Lagrangian Eulerian (ALE) coordinates in fluid mechanics. 

The moving domain 17(f) will be represented as the image of a time-dependent family of diffeomorphisms T(/:): 
17 1 —> 17(f). Let N represent the outward pointing unit normal to V and let T(f) be given by 

T(t) = {x\ x = xo + h(t,Xo)N, x 0 er}. 

Assuming that the signed height function h(t,-) is sufficiently regular and F(f) remains a small graph over F, we can 
define a diffeomorphism : 17 —► 17(f) as the elliptic extension of the boundary diffeomorphism xq xq + h(t,xo)N, 
by solving the following Dirichlet problem: 

AT = 0 in 17, (1.4) 

H7(f,x) =x + h(t,x)N(x) x&Y. 

We introduce the following new variables set on the fixed domain 17: 


q=p otp 


v = — Vpo vp 
A= [Z7\P] -1 
J = deti7T' 


{temperature), 

{“velocity”), 

{inverse of the deformation tensor), 
{Jacobian determinant), 


We now pull-back the Stefan problem (ll.lt from 17(f) onto the fixed domain 17. If we let g denote the Jacobian of 
the transformation H7(f,-)|r : F — >T(t), and let n(t,-) denote the outward-pointing unit normal vector to the moving 
surface F(f), then the following relationship holds fl5l : 

1 \fg Tli o vp(f, x) = A\ (f , x)N k (x). 

It thus follows that the outward-pointing unit normal vector n(t,-) to the moving surface F(f) can be written as 
(no'fy)(t,x) = A T N/\A T N\ . We shall henceforth drop the explicit composition with the diffeomorphism 4>, and 
simply write 

n(t,x) = A T N/\A T N\ 

for the unit normal to the moving boundary at the point tp( t,x ) € T(f). 
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The classical Stefan problem on the fixed domain fl is written as (see 112711281 ) 


& 

II 

•Os 

in [0 

(1.5a) 

O 

II 

in [0,T) x f2, 

(1.5b) 

q = 0 

on [0,T) x T, 

(1.5c) 

u v-A t N 

* — N-A T N 

on [0,T) x T, 

(1.5d) 

ATi = 0 

on [0,T) x fl, 

(1.5e) 

vp = e + hN 

on [0,T) x T, 

(1.5f) 

q = q o >0 

on {t = 0} x 12, 

d.5g) 

h = 0 

on {t = 0} x T, 

(1.5h) 


Problem (11.51 ) is a reformulation of the problem (11.11 ). Observe that the boundary condition (1 1 .5db is equivalent to 


' i t -n(t)=v-n(t ) on [0,T)xT so that \I/(f)(r) =T(f), (1.6) 

which is but a restatement of the Stefan condition ( ll.lbb . Since the factor N- A T N will show up repeatedly in various 
calculations, it is useful to introduce the abbreviation: 

A = N-A t N. (1.7) 


Note that initially A = 1 and it will remain close to 1, since for small h the transition matrix A remains close to the 
identity matrix. 

Since the identity map e : O —> Q is harmonic in Q and T — e = hN on F, standard elliptic regularity theory for 
solutions to (11.41 ) shows that for t £ [0,T), 

II^V) - e|| ff s(n) <C l ||/i(f,-)|| ffs -°.5( r ), s>0.5, 

so that for h sufficiently small and s large enough, the Sobolev embedding theorem shows that VT is close to Id, and 
by the inverse function theorem, 4/ is a diffeomorphism. 


1.3.1. The high-order energy and the high-order norm. We will specialize to the case <7 = 2 for the remainder of this 
paper. The case d = 3 requires only our norms to contain one more degree of differentiability, while the rest of the 
argument is entirely analogous. 

To define the natural energies associated with the main problem, we must employ tangential derivatives in a neigh¬ 
borhood which is sufficiently close to the boundary F. Near l = di 2, it is convenient to use tangential derivatives 
d a , while away from the boundary, Cartesian partial derivatives d Xi are natural. For this reason, we introduce a 
non-negative C°° cut-off function //: O —> R_ with the property 


fi(x) = 0 if \x\<p] p{x)=l if dist(a;,r) <cr. 

Here p,a SK + are chosen in such a way that B p ( 0) <s tt and (a;|dist(a:,r) < cr} £ fl\f? p (0). 

Definition 1.1 (Higher-order energies). The following high-order energy and dissipation functionals are fundamental 
to our analysis: 


£{t) = £{q,h){t) = 

\ E \\T 1/2 d S d b t v\\h x + 1 - 

|S|+2b<5 

E \\{i-p) l/ 2 d s d b t v\\li 

|<S|+2h<5 


E + i E \\T 1/2 (d s d»q + d s d b t*- v )\\ 

|5|+2b<6 \a\ + 2b<6 

+ ^ E m-ri 1 / 2 (d 3 d b t q + d s d^-v)f Ll 

|5|+2b<6 


2 

L l 
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and 

V(t)=V(q,h)(t)^ 

E IIm 1/2 ^M||,+ E \{-d N q) 1 / 2 Ad s d^ t \l 2x + E llM 1/a (^ft+^#®t-t/)||£, 

|a|+2b<6 |5|+2fc<5 |a|+26<5 

+ e ii(i-m) 1/2 ^mh 1+ e ii(i-M) 1/a (^ft+fl a ^f«)ii^, 

|a|+26<6 |a|+2b<5 

where we recall the definition of A given in (11.7b . Finally, we introduce the total energy E(t): 

E(t) = sup £(t) + f V(r)dT. (1.8) 

0<s<£ Jo 

Note that the boundary norms of the gauge function are weighted by -f—d^q. We thus introduce the time- 
dependent function 

X(t) = inf (-d N q){t,x)>0, 

xGF 

which will be used to track the weighted behavior of h. It is important to note, that due to the smoothness assumption 
on r it is easy to see that for any local coordinate chart (c) Sl ,..., d Sd _ 1 ) for F we have the equivalence 

E |(-3 JV g) 1/2 A0 a 0**|£ s « E \(-dNq) 1 / 2 d^...d^h\ 2 L2{r) , (1.9) 

|S|+2fc<6 ii=03i.-,/3d_i) 

|/3|+2b<6 

where X « Y means that there exist positive constants C\ and C 2 such that C\Y < X < CfY. In our case, the two 
constants depend on the choice of the local chart. 

Definition 1.2 (High-order norm). The following high-order norm is fundamental to our analysis: 

3 2 

— E + WqWl^H 6 - 5 +y^Mdt. ( lt\\ 2 L 2 H 5-2l 

1=0 1=0 

+ sup e l3s \\q{sc)\\ 2 H 5+ E 

°- s - 4 |5|+2i<6 

3 2 

+ x(t)J 2 \d l t h\ 2 LoaH 6 . 2 l +x{t)j 2 \d l t + 1 h\ 2 L 2 H 5-2l + \h\^aa H 4.5 
1=0 1=0 

(1.10) 

Flere /3 = 2A — 77 , where A is the smallest eigenvalue of the Dirichlet-Laplacian on f l and 17 >0 is a small but fixed 
number to be determined later. 

Remark 1.3. A subtle feature of the above definition is the loss of a ^-derivative-phenomenon for the temperature q. 
By the parabolic scaling (where one time derivative scales like two spatial derivatives), one might expect q to belong 
to L 2 iF 7 ([0,T);H), since d l t + 1 q£L 2 H 5 ~ 2 l ([0,T)]fl), for 1 = 0,1,2. This is, however, not the case, as the height- 
evolution equation ( 11.5db scales in a hyperbolic fashion, and thus places a restriction on the top-order regularity of 
the unknown q , allowing only for q£L 2 H 65 ([0,T);Q). 

1.4. Steady states. Note that any C 1 simply connected domain represents a steady state of (11.1b . In other words, 
for any simply connected domain OsC 1 , the pair (u = 0,f = Oil) forms a time-independent solution to (11.1b . In 
particular, it is challenging to determine which steady state a small perturbation will decay to. Thus the problem 
of asymptotic stability, rather than the optimal regularity of weak/viscosity solutions, is one of the main motivating 
questions for this work. In particular, we work with classical solutions with a high degree of differentiability on the 
initial data. 

1.5. Rayleigh-Taylor sign condition or non-degeneracy condition on q Q . With respect to qo, condition (fO > be- 

inf [—<9 at<7o (a;)] > 5 > 0 on T. 
a;Gr 


comes 
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For initial temperature distributions that are not necessarily strictly positive in f l, this condition was shown to be 
sufficient for local well-posedness for (fO (see ll27l [3911411 ). On the other hand, if we require strict positivity of our 
initial temperature functiorfl, 

qo>0 infl, (1-11) 

then the parabolic Hopf lemma (see, for example, j20l ) guarantees that —d^qU, x) > 0 for 0 < t < T on some a priori 
(possibly small) time interval, which, in turn, shows that £ and V are norms for t > 0, but uniformity may be lost as 
t —>-0. To ensure a uniform lower-bound for —diyq(t) as t—y 0, we impose the Rayleigh-Taylor sign condition with 
the following lower-bound: 

— dNqo > C* / qotpidx, (1.12) 

J n 

Here, ip\ is the positive first eigenfunction of the Dirichlet Laplacian —A on f2, and (7* >0 denotes a universal 
constant. The uniform lower-bound in (11.12b thus ensures that our solutions are continuous in time; moreover, (11.12b 
allows us to establish a time-dependent optimal lower-bound for the quantity \(t) = infxer(—dNq){t,x) > 0 for all 
time t > 0, which is crucial for our analysis. 


1.6. Main result. Our main result is a global-in-time stability theorem for solutions of the classical Stefan problem 
for surfaces which are assumed to be close to a given sufficiently smooth domain O and for temperature fields close to 
zero. The notions of near and close are measured by our energy norms as well as the dimensionless quantity 


11 gQ 11 4 

IMIo' 


(1-13) 


as expressed in the following 


Theorem 1.4. Let ( qo,ho ) satisfy the Rayleigh-Taylor sign condition ( 17.721) . the strict positivity assumption ( 17.771) . 
and suitable compatibility conditions. Let K be defined as in ( I7.7ib . Then there exists an cq > 0 and a monotonically 
increasing function F : (l,oo) —> R+, such that if 


5(0) < 


£Q 

F(i<y 


then there exist unique solutions ( q,h ) to problem ( 17.51) satisfying 


(1.14) 


S(t)<Ce o, f€[0,oo), 

for some universal constant C > 0. Moreover, the temperature q(t) —> 0 as t —> oo with bound 

where /3 = 2X — O(eo) and A is the smallest eigenvalue of the Dirichlet-Laplacian on LI. The moving boundary T(t) 
settles asymptotically to some nearby steady surface T and we have the uniform-in-time estimate 

sup |h(t,-)-/i 0 |4.5<\/eo 
0<£<oo 


Remark 1.5. The increasing function F(K) given in ( 17.74b has an explicit form. For universal constants (7,(7 > 1 
chosen in Section [7] 

F(K) = matx{8K 2ceK \c 10 (\nK) w K 2 ° ex }. (1.15) 

Remark 1.6. The use of the constant K in our smallness assumption ( 17.741 ) allows us to determine a time T = I f 
when the dynamics of the Stefan problem become strongly dominated by the projection of q onto the first eigenfunction 
ipi of the Dirichlet-Laplacian. Explicit knowledge of the K-dependence in the smallness assumption ( 17.74b permits 
the use of energy estimates to show that solutions exist in our energy space on the time-interval [0 ,Tk\- For t > Tk, 
certain error terms (that cannot be controlled by our norms for large t) become sign-definite with a good sign. 

Remark 1.7. An analogous theorem was stated in 11281 . for perturbations of steady surfaces initially close to a sphere. 
Therefore, this work generalizes that result. Moreover, our methods are general enough to apply to other geometries 
as well. An example is that of a free boundary parametrized as a graph over a periodic flat interface. 


'Condition Gil} is natural, since it determines the phase: Q(t) = {q(t) > 0}. 
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Remark 1.8 (On compatibility conditions). The first compatibility condition on the initial temperature q o is 

9o|r = 0. 

The second condition arises by restricting the parabolic equation (1 1.5 ab to the boundary T and using the boundary 
conditions (1 1 .5cb and (11.6b . It gives 

dNN<lo + {d~- l)n r dNq 0 + (d]s[qo) 2 = 0 on T. 

Here Kr stands for the mean curvature ofT. Higher order compatibility conditions arise by taking time derivatives 
of d 1 .5ab , re-expressing them in terms of purely spatial derivatives via (11.5ab and restricting the resulting equation to 
the boundary T at time t = 0. 

Remark 1.9. An interesting problem is to determine the asymptotic attractor - the steady state 1 just from the initial 
data (u 0 ,r 0 ). This is strongly connected to the so-called momentum problem, which is a problem of determining the 
domain Qfrom the knowledge of its harmonic momenta c^ = (fr.M. d — >R, A^> = 0. A related question arises 

in the Hele-Shaw problem, see 0261 . 

1.7. Local well-posedness theories. In |[27l , we established the local-in-time existence, uniqueness, and regularity 
for the classical Stefan problem in /7-bascd Sobolev spaces, without derivative loss, using the functional framework 
given by Dehnition ll.il This framework is natural, and relies on the geometric control of the free-boundary, analogous 
to that used in the analysis of the free-boundary incompressible Euler equations in 1(1411151 ; the second-fundamental 
form is controlled by a a natural coercive quadratic form, generated from the inner-product of the tangential derivative 
of the cofactor matrix JA, and the tangential derivative of the velocity of the moving boundary, and yields control of 
the norm f r (—d]yq(t))\d k h\ 2 dx' for any k > 3. The Hopf lemma ensures positivity of —d]\rq(t) and the Taylor sign 
condition on qo ensures a uniform lower-bound as t —> 0. 

The first local existence results of classical solutions for the classical Stefan problem were established by Meir- 
manov (see l(39l and references therein) and Hanzawa |[29l . Meirmanov regularized the problem by adding artificial 
viscosity to (ll.lbb and fixed the moving domain by switching to the so-called von Mises variables, obtaining solutions 
with less Sobolev-regularity than the initial data. Similarly, Hanzawa used Nash-Moser iteration to construct a local- 
in-time solution, but again, with derivative loss. A local-in-time existence result for the one-phase multi-dimensional 
Stefan problem was proved in Il24l , using /A-type Sobolev spaces. For the two-phase Stefan problem, a local-in-time 
existence result for classical solutions was established in ED in the framework of //'-maximal regularity theory. 

1.8. Prior work. There is a large amount of literature on the classical one-phase Stefan problem. For an overview we 
refer the reader to lf22l(39l |46 1 as well as the introduction to || 28 1 . First, weak solutions were defined in 0111211. [37 1. For 
the one-phase problem studied herein, a variational formulation was introduced in 11231 . wherein additional regularity 
results for the free surface were obtained. In (6) it was shown that in some space-time neighborhood of points xq 
on the free-boundary that have Lebesgue density, the boundary is C 1 in both space and time, and second derivatives 
of temperature are continuous up to the boundary. Under some regularity assumptions on the temperature, Lipschitz 
regularity of the free boundary was shown in 0- In related works mm it was shown that the free boundary is 
analytic in space and of second Gevrey class in time, under the a priori assumption that the free boundary is C 1 with 
certain assumptions on the temperature function. In J9) the continuity of the temperature was proved in d dimensions. 
As for the two-phase classical Stefan problem, the continuity of the temperature in d dimensions for weak solutions 
was shown in Eo). 

Since the Stefan problem satisfies a maximum principle, its analysis is ideally suited to another type of weak solu¬ 
tion called the viscosity solution. Regularity of viscosity solutions for the two-phase Stefan problem was established 
in a series of seminal papers 00. Existence of viscosity solutions for the one-phase problem was established in (32), 
and for the two-phase problem in lf33) . A local-in-time regularity result was established in mu, where it was shown 
that initially Lipschitz free-boundaries become C 1 over a possibly smaller spatial region. For an exhaustive overview 
and introduction to the regularity theory of viscosity solutions we refer the reader to El. In |36) the author showed 
by the use of von Mises variables and harmonic analysis, that an priori C 1 free-boundary in the two-phase problem 
becomes smooth. 

In order to understand the asymptotic behavior of the classical Stefan problem on external domains, in 142) the 
authors proved that on a complement of a given bounded domain G, with non-zero boundary conditions on the fixed 
boundary dG, the solution to the classical Stefan problem converges, in a suitable sense, to the corresponding solu¬ 
tion of the Hele-Shaw problem and sharp global-in-time expansion rates for the expanding liquid blob are obtained. 
Moreover, the blob asymptotically has the geometry of a ball. Note that the non-zero boundary conditions act as an 
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effective forcing which is absent from our problem and the techniques of l42l do not directly apply. Since the cor¬ 
responding Hele-Shaw problem (in the absence of surface tension and forcing) is not a dynamic problem, possessing 
only time-independent solutions, we are not able to use the Hele-Shaw solution as a comparison problem for our 
problem. 

A global stability result for the two-phase classical Stefan problem in a smooth functional framework was also 
established in |39l for a specific (and somewhat restrictive) perturbation of a flat interface, wherein the initial geometry 
is a strip with imposed Dirichlet temperature conditions on the fixed top and bottom boundaries, allowing for only 
one equilibrium solution. A global existence result for smooth solutions was given in ED under the log-concavity 
assumption on the initial temperature function, which in light of the level-set reformulation of the Stefan problem, 
requires convexity of the initial domain (a property that is preserved by the dynamics). 

Remark 1.10. We remark that global stability of solutions in the presence of surface tension does not require the 
use of function framework with a decaying weight, such as —di\/q(t). In this regard, the surface tension problem is 
simpler for two important reasons: first, the surface tension contributes a positive-definite energy-contribution that 
is uniform-in-time, and provides better regularity of the free-boundary (by one spatial derivative), and second, the 
space of equilibria is finite-dimensional and thus it is easier to understand the degrees-of-freedom that determine the 
asymptotic state of the system. 

1.9. Methodology. Broadly speaking, our methods combine high-order energy estimates with maximum principle 
techniques. Once the problem is formulated on the fixed domain with the help of the harmonic gauge explained above, 
we notice that the natural quadratic energy quantities that track the regularity behavior of the moving boundary, come 
weighted with the normal derivative of the temperature. This weight is a time-dependent quantity and its evolution 
is tied to the free boundary itself. This coupling is nonlinear and it is one of the central difficulties in closing our 
estimates. 

Our strategy is based on ll28l and it contains three basic steps. We first show that under the assumption of smallness 
on the norm S(t) over some time interval [0,T], the energy E and the norm S are equivalent, i.e. 

S(t)<E(t)<S(t), te [0,T], (1.16) 

Our second step is to establish the key energy inequality in the form 

E(t)<C 0 +i ^ [ (d N q t )\d s d l t h\ 2 dS(T) + P(S(t)), (1.17) 

|a|+2Z<6 1 ' 0 

where P is a cubic polynomial (see (11.3b ) and Co is a small quantity depending only on the initial data. Combin¬ 
ing (11.16l > and (11.17b . we infer that 

S(t)<C 0 + C J2 I f ( d N q t )\d 5 d l t h\ 2 dS(T)+P(S(t )) (1.18) 

|a|+2Z<6^° 

'-v-' 

dangerous term 

on the time interval of existence. If it were not for the sum on the right-hand side above, a simple continuity argument 
would yield a global existence result for small initial data. However, the sum appearing on the right-hand side of (11.1 8b . 
while seemingly cubic, cannot be bounded by P(S(t)). Instead, in the third step we show that after a certain, precisely 
quantified amount of time, this “dangerous term” becomes negative and can thus be trivially bounded from above by 
zero. 

The key novelty with respect to If28ll is a new quantitative lower bound on the weight —d^q which appears in our 

definition of the energy E(t). Note that this quantity is expected to converge exponentially fast to 0 as the unknowns 

settle to an asymptotic equilibrium. We employ the theory of “halfeigenvalues” associated with the Bellman-Pucci- 
type operators to generate a comparison function, which then allows us to use the maximum principle and get a nearly 
sharp lower bound: 

-< 9 Ar<?>e(- A+ °( £))t , 

where A denotes the first Dirichlet eigenvalue associated with the domain O. In our previous work H28| , we relied on 
a rather explicit Bessel-type comparison functions used by Oddson in flOl . which in particular, required that we work 
in a nearly spherical domain. The above lower bound is much more flexible and it is explained carefully in Section [3] 

The presentation in the paper is considerably simplified with respect to ll28l and we believe that our energy method 
in conjunction with maximum principles can be useful for the stability analysis in other free boundary problems in 
absence of surface tension. 
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1.10. Plan of the paper. In Section [2] we introduce the bootstrap assumptions and formulate the equivalence rela¬ 
tionship between the energy and the norm. In Section [3] we provide a dynamic lower bound estimate on \(f). This is 
the main new ingredient with respect to |f28l and we use the theory of half-eigenvalues for the Pucci operators. Finally, 
in Section Q] we give the proof of Theorem 1 1.41 thereby explaining our continuity method as well as a comparison 
argument used to show the sign-definiteness of the “dangerous linear terms” described above. 


2. Bootstrap assumptions and norm-energy equivalence 

2.1. The bootstrap assumptions. Let [0,T) be a given time-interval of existence of solutions to (1 1 ,5b . We assume 
that the following two assumptions hold: 

S(t) < e, i€[0,r), (2.19) 

X(t)>c t€[0,T), (2.20) 

where e and r/ are to be chosen sufficiently small later and A stands for the first Dirichlet eigenvalue associated with 
the domain f 2 . 

2.2. Norm S and total energy E are equivalent. Recall the notation introduced in ( 11.91 ). 

Proposition 2.1. There exists a sufficiently small e 1 such that if S{t) < e' on a time interval [0,T] then 

S(t)~E(t), Vt€[0,T]. 

Proof. The proof of this fact is one of the pillars of our strategy. It has been presented in detail in Sections 2.1 - 2.5 
and Section 4.2 of ll28l and, therefore, we omit it here. We note that the direction 5(f) < Eft ) is obviously harder to 
prove, as the energy function Eft) a-priori controls only tangential derivatives of the temperature q. In ll28l we use 
a version of the elliptic regularity statement for equations with Sobolev-class coefficients to obtain control of normal 

derivatives (see ED). □ 


3. Lower bound on x(i) and improvement of the second bootstrap assumption 
The heat equation ( 1 1 .5ab for q can be written in non-divergence form as 

(3.21a) 
(3.21b) 
(3.21c) 


qt-akjq,kj~hq,k = o in 12, 

<7 = 0 on r, 

< 7 ( 0 ,•) = <70 > 0 in 12 


where the coefficient matrix a = (akj)k,j= i,2> and the vector b = (for, bf) are explicitly given by 


O-kj 




b^A^Al+Af^. 


(3.22) 


By the bootstrap assumption (12.191 ) and the definition ( 11.101 ) of 5(f), we have that \h\i .5 < yfi on [0,T), and there¬ 
fore by the Sobolev embedding H 1 ( r) °->-L 00 (r), we infer that \h\w^.°° \T ( - From this observation, (13.22b . and the 

definition of the transition matrix A, we infer that 


\ a kj Skj \ fS \/e, — 1)2), 

\bi\<V~e, {i = 1,2). 


Therefore, there exists a constant K > 0 such that the ellipticity constants associated with the matrix (oij)* j=i,2 are 
between the values fi[ = 1 — and fi ' 2 = 1 + uniformly over [0,T). 

Before we proceed with calculating a lower bound for x(f), we briefly explain the Bellman operators ||2l [5l [l9l 
l25l(38l which are closely connected to the well-known extremal Pucci operators. They will allow us to formulate a 
nonlinear analogue of the “first” eigenvalue for the elliptic part of the operator defined in (13.21al) . 

Let 12 be an arbitrary simply connected C 1 -domain. We define the extremal Pucci operator ; , 2 l[25l [5l with 
parameters 0 < /ri < 7 x 2 by 

= rJnf C T( X )- 


( 3 . 23 ) 
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Here /C MllAt2 denotes the set of all linear second-order elliptic operators, whose ellipticity constant is between Mi and 
M 2 , i-e., 

A^ ilj ^ 2 ={L| L = a-ijdij +bidi + c, (3.24) 

feR d }. 

It is well known that the operators Mf, /l2 are, in general, fully nonlinear second-order elliptic operators, positive, 
and homogenous of order one. The latter property allows us to formulate an associated “eigenvalue” problem, looking 
for the solutions of 

~M~ 1 ^ 2 u = Xii in n, (3.25) 

u = 0 on dfl. 

We next state some of the results from f38l that that will play an important role in this paper (for further references on 
the so-called half-eigenvalues associated with positive homogenous fully nonlinear operators we refer the reader, for 
example, to |[5l[2lH9ll): 

• There exist two positive constants Ai and A 2 called the first half-eigenvalues and two functions o-\. 0-2 4 
C 2 (fi)nC'(fi) such that (Ai,pi) and (X 2 ,g 2 ) solve (13.251) . and £>1 > 0, g 2 < 0 in il. 

• The first two half-eigenvalues are simple, i.e. all positive solutions to (13.25b are of the form (Ai,api) with 
a > 0 and analogously, all negative solutions are of the form (X 2 ,aQ 2 ), « > 0. 

• Finally, the first two half-eigenvalues are characterized in the following manner: 

Ai= sup p{A), X 2 = inf p(A), (3.26) 

AgK i _, 1 ,^, 2 AeK .^ 1,^2 

where // ( A) stands for the smallest Dirichlet eigenvalue associated with the second order linear elliptic oper¬ 
ator A. 

3.1. Lower bound on x(f) and the improvement of (12.201 ). The key ingredient to the proofs of Propositions 12.11 
and 14. II is a quantitative lower bound on the weight x(t)- This is achieved by using the maximum principle and 
constructing an appropriate comparison function. 

Lemma 3.1. Under the bootstrap assumptions ( 12.191 ) - (12.20b with e sufficiently small, the following inequality holds: 

where c\ = f^qoifiidx is the first coefficient in the eigenfunction expansion of the initial datum go with respect to the 

L 2 orthonormal basis {<Mi, <M 2 , ■ • •} of the eigenvectors of the operator — A on f1, i.e go = C\p\ -\-c 2 P 2 + Moreover, 

A stands for the smallest Dirichlet eigenvalue associated with the domain fl and A (t) satisfies the estimate: 

|A(t)| <C\fe. 

In particular, with e > 0 sufficiently small so that Cy/e < rj/4, we obtain the improvement of the bootstrap bound (12.20b 
given by x{t) > cie - ( Al+?? / 4 ) t . 

Proof. Let us choose mi = 1 — Kyfe and p, 2 = l + Ky/e. Recall that K was defined in the paragraph after (13.22b . It 
follows that L € /C ; ,, >M2 . We let p i be the first half-eigenvector associated to /l2 as above. Consider the following 
comparison function 

v(t,x) = e q 1 . 

Note that v vanishes on Oil = T. A straightforward calculation together with the definition of ,M ( , jL2 shows that 

( d t — L)v = — X\V — e -Al t Lgi 

<-Aiu-e“ Alt M" i/i2 pi 
— —Aiu + e Al< Aipi 
= 0 . 

Therefore v is a subsolution to the parabolic problem ( 13.21b . The next key observation is that the eigenfunction pi (a;) 
behaves like a constant multiple of the distance function dist(x, Fj as x approaches the boundary T. Namely, since the 
operator is concave, the solution is C 2, “ 11441 (81 and the Hopf lemma —d\-g\ > 0 holds (see for instance Lemma 
2.1 in 0). Therefore, function v behaves like c dist(.x-,L)e A i 4 as x approaches the boundary T for some constant c. 
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Here dist(x,F) denotes the distance function to the boundary I’. We first want to show that for any arbitrarily small 
time a > 0 there exists a strictly positive constant 5(a) > 0 such that q — 5v is a positive supersolution to the parabolic 
problem ( 13.211 ) on the time interval [a,T). 

Since v is a subsolution and q is a solution, it follows that for any 5 > 0, q — 5v is a supersolution. The positivity of 
q — 5v at t = a follows from the parabolic Hopf lemma, from which we infer the existence of a constant 5(a) such that 
^ > 5(a) uniformly over Cl. Note that we have used the fact that v(a,x) behaves like c x dist(x) near the boundary T 
for some positive constant c. Thus by the maximum principle, q — 5(a)v> 0 on [a,T). This implies 

q(t,x) > 6(a)v(t,x) > CJ(cr)dist(a;,r)e _Alt , t £ [a,T), 


which yields 


dq(t,x) 

dN 


>C5(a)e~ Xlt , te[a,T). 


The above estimate is however not yet satisfactory, as the constant 5(a) may degenerate as cr goes to zero. 
We now revisit our usage of the parabolic Hopf lemma above. For small t > 0 let 

O t = {xeO| dist(x,T)>f}, t> 0. 


Note that f l t is a compact proper subset of Cl. From the proof of the parabolic Hopf lemma (see for instance Theorem 
3.14 in l20l ). the value —dq/dN | t=(T is proportional to the minimal value of the temperature q on a space-time region 
strictly contained in the space-time slab K t :=Cl t x [f/2,37/2] x [0,2t] divided by t (which is proportional to the 
distance of K t from the parabolic boundary of 1 1 x [0,2t]). Note that, as t approaches 0 we may loose uniformity-in- 
time in our constants. This is however not the case since <9y q is continuous at t = 0 and by the assumption (11.121 ) 

-d N qo= 9n<10 ci >C'*Ci. (3.27) 

Cl 

Assumption ( 11.121) is used only in (13.271 ) to insure that there exists a universal constant C* independent of ci such that 
L = ( iTvf/o)/ci > C t . The quantity L is dimensionless, and the assumption L > C* is not a restriction on the initial 
data. In other words, if we had not assumed (1 1.1 2b . the only modification in the statement of the main theorem would 
be that the smallness assumption on initial data (11.14b is additionally expressed in terms of L as well. 

As to the bound on A, note that by (13.26b . the exponent Ai is characterized by the condition 


Ai= sup p(A). 

Zt£^C/Al,/i 2 


Since |£ij —1| <-\/e, * = 1,2, it follows that for any matrix Ae/C Mli/i2 the estimate |— Id| < -^/e holds. Since the 
function p(-) is a continuous function from the space of 2 x 2 matrices into R, it thus follows that 

|A| = |Ai-/z(Id)| = | sup p(A)-p(ld)\<y/e. 


4. Energy estimates and improvement of the first bootstrap assumption 
Proposition 4.1. Assuming the bootstrap assumption (12.19b and with e > 0 chosen sufficiently small, 

E(t)<Co +i ^ r [(-d N q t )\d s d l t y\ 2 dS(r) + CP(S(t)), (4.28) 

|a|+2Z<6"'° 

where Cq depends only on the initial data, C >0 is a generic positive constant depending only on the dimension d, 
and P denotes an order-r polynomial with r > 3 of the form (Ob . 

Proof. The proof of the proposition is entirely analogous to the proof of Proposition 3.4 from 128). □ 

Proposition 4.2. Let the solution (q,h) to the Stefan problem ( 17.5b exist on a given maximal interi’al of existence 
[0,7”) on which the bootstrap assumptions ( 12.79b and (12 .20b are satisfied. 

(a) There exists a universal constant C such that if the smallness assumption ( 17.74b for the initial data holds and 
ifT> T k = C In K, then 

- qt(T K ,x)>Cc 1 e- X ' TK 


pi(x), xeB i(0) 
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where ipi is the first eigenfunction of the Dirichlet-Laplacian on Q and c\ = J^qopi dx. As a consequence, 

mfd N q t (T K ,x)>0. 

(b) With the smallness assumption ( 11.141 ). we indeed have the bound T > C\wK. 

(c) Moreover, under the same assumption as in part (b), the following lower bound on d^q(t.x) holds: 

inid N q{t,x)>Q, t£[T K ,T). (4.29) 

Proof. The proof of part (a) of is the same as the proof of Lemma 4.2 in ll27l . 

As to the proof of part (b), we start by making the claim that the dangerous term from the inequality ( 14.281 ) satisfies 
the bound 


|5|+2i<6 

Note, that if |a| + 2 1< 6, then 
/■* 

(— On qt ) I d a d\ T 1 1 z dS dr 


V f [(-d N q t )\d s d l t 'S>\ 2 dS{F) <CK 2 f e I ' T S'(r)dr. 

, o Jr J o 


(4.30) 


to 


In order to bound the term 


dNqt (-d N q) \d s d l M 2 


dSdr<C 


dwqt 


dNq 


S(r)dT. 


—dNq 

| |, we need a decay estimate for the numerator |djv<Zt|- The Sobolev embedding theory 

would yield the bound \dNqt\ l°° < ||< 7 t|| 2+<5 for S > 0, but by definition of the norm S, it is only the fJ 2 (fl)-norm of q t 
for which we have the desired decay. We obtain the decay estimate for q t from Appendix B of | 

\d N qt\L^<K 2 Cl e-^ 2 . 

It then follows from the bootstrap assumption (12.201 ) that 

CK 2 Ci e _ ( Al_,? /2) T 
° 1 < CK 2 e r,T . 


(4.31) 


dNqt{r) 


d N q(r) 


< 


c 1 e-( Xl+r >/ 2 ') T 


which, in turn, establishes (14.301) . In conjunction with Proposition ^. 1 1 this yields the bound 


E{t)<E(0) + CK 2 f e nT S( t) dr+ CeS(t). 
Jo 


By Proposition ^. II with e sufficiently small, we conclude that 


E(t)<2E{0) + CK 2 [ e VT E{ t) dr, t€[0,T\, 
Jo 


(4.32) 


(4.33) 


where T is the maximal interval of existence on which the bootstrap assumptions (12.19b and (12.20b hold (with e 
sufficiently small). A straightforward Gronwall-type argument based on (14.33b . identical to Step 1 of the proof of 
Theorem 1.2 in l27l . implies that as long as the rj from the bootstrap assumption (12.20b is smaller than 6'In K, the 
maximal interval of existence [0,T), on which both the bootstrap assumptions (12.19b and (12.20b are valid, satisfies 
T > ChiK, and the following exponentially growing bound holds: 

E{t)<2E(f))e CK2 \ t£ [0,7"). (4.34) 

To prove the part (c), we resort to maximum principle techniques once again. To this end, we define a barrier 
function if> to be the solution of the following elliptic problem 

A%fj=—1 in ft (4.35) 

ip = 0 on T. 

We then define the comparison function T : |0,T) via 

E(t,x) = Kie~^ xt (ipi(x) — K 2 i/t), (4.36) 


with positive constants ki,K 2 to be specified later. A straightforward calculation shows that 

3 1 3 

(d t - atjdij - bidf)T = K X e~ * xt [- -\<pi -k 2 + - {atj - 8ij)(ip i - n 2 f>) -h-(Vp i - n 2 V^)]. (4.37) 

Note that the first and the second terms in the parenthesis on the right-hand side of ( 14.37b are negative, while the fourth 
and the fifth terms, are small, being of order e. If x is close to T, then the second term dominates the third term and if 
x is away from the boundary T, then one can choose n 2 > 0 so that the first term dominates the third term. Thereby we 
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use the fact that ip i and ib both vanish at F, they are both non-negative (by the maximum principle), and both satisfy 
the Hopf lemma (since they are both super-solutions). It follows, then, that there exists a /t 2 > 0 and some constant C\ 
such that 

(dt — dijdij — bidi)bF < —CiKie~^ xt . (4.38) 

It then follows from (14.381) and (13.221) that 

(d t - dijdij -bidi)(-q t - F) > -(dtaijq,ij+d t b i q i + d t A k z q, k w l + A/q, k w l t ) + C 1 n 1 e~? xt . (4.39) 

Note, however, that the term in parenthesis on the right-hand side above is a quadratic non-linearity and as such decays 
at least as fast as e - ^*: 

\\dtCHjqnj+dAqi + dtA^q^w 1 + A^q, k w l t \\ L ^ <C 2 c 1 ee~^ t . 

Now, using ( 14.39b and the above bound, we note that by choosing the constant k\ = pfcq e, we have that 

(d t - dijdij - bidi)(-q t -bF)> C 2 ciee~i xt - C 2 c\ee~ pt > 0, 

since /3 = 2 A — 77 > | A. The previous bound implies that — q t — T is a supersolution for the operator dt. — a,;, d l3 —bidi. 
Moreover, by the construction of bF , we have -qt — F 0 on F. Furthermore, at time Tk = 6”hi K, we have by the 
part (b) of the proposition and (14.361) . that 

(-qt - bF)\ T =c lnK > Cc\e~ XT ipi(x) - Cciee^^ XT ipi ( x ) + Caen 2 e~i XT ^(x) > 0 

for e sufficiently small. Thus, as in the proof of Lemma l3TTl there exists a constant m > 0 such that 

—q t (t,x) — F(t,x) > m dist(a;,r)e _ ^ A+0 ^^ t , t > Tk, 


or, in other words, 


ip(x) 


-„((,*) > mdi s t(x,r) e -(»°(-»* + Cc ie dist(x,r)e (5^) -^ dist(l>r) 

= (m + Cci «<-*«-<>«))• (£gh - ). 

which readily gives the positivity of dwqt on the time-interval [Tk,T[ since di^fr7 > 0 by our choice of 

k 2 above. We conclude that the positivity of —q t at time Tk = C'lu K is a property preserved by our bootstrap regime 
and, moreover, we obtain a quantitative lower bound on d^qt on the time interval \Tk,T\. □ 


Remark 4.3. In the proof of part (b) of Proposition \4~2\ we made a rather crude use of the energy estimate given by 
Proposition 1 4. 1\ In particular, we cannot use this argument to prove global existence, as the constants grow in time; 
however, in part (c) of the proposition, we have used a more sophisticated argument based on the maximum principle 
to infer the sign-definiteness of the term d^qt after a fixed amount of time has passed. 


Proof of Theorem 1 1.41 Assume for contradiction that T < oo. For any t £ the energy identity takes the form 

E (t)+^Y,j T j r dNqt\d 3 d l t ^\ 2 dS<E(T K ) + P(S(t))<E(T K ) + 0(e)E(t). 

Note here the absence of the exponentially growing term in the above bound as compared to the inequality (14.34b . This 
is due to the fact that terms J T ^ f r djvqt\d a d l t ^\ 2 dx, |5| +21 < 6, are positive and no longer treated as error terms. By 
absorbing the small multiple of E(t) into the left-hand side, and using the positivity of dsj q f from Step 2, we obtain 
that 

E(t)<2E(T K )<8E(0)e 2CK2TK , t£[T K ,T), (4.40) 

by (14.34b . Finally, we choose eo in the statement of Theorem ll.4l so that eo < e/2. The bound ( 14.40b and the condition 
E( 0) < e o/F(K) (with F(K) given as in (11.15b ) imply 

t€[T K ,T). 

Together with Lemma I3TT1 we infer that the bootstrap assumptions ( 12.19b and (12.20b are improved. Since E(-) is 
continuous in time, we can extend the solution by the local well-posedness theory to an interval [0,T +T*] for some 
small positive time T*. This however contradicts the maximality of T if T were finite and hence T = oc. This 
concludes the proof of the main theorem. 




























GLOBAL STABILITY OF STEADY STATES IN THE CLASSICAL STEFAN PROBLEM 


13 


Acknowledgements 

SS was supported by the National Science Foundation under grant DMS-1301380 and by the Royal Society Wolfson 
Merit Award. MH was partly supported by the National Science Foundation under grant DMS-1211517. Some of this 
work was completed during the program Free Boundary Problems and Related Topics at the Isaac Newton Institute for 
Mathematical Sciences at Cambridge, UK. We are grateful to the organizers, Gui-Qiang Chen, Henrik Shahgholian 
and Juan Luis Vazquez, for both the invitation to participate in the program and to contribute to this special volume. 

References 


[1] Adams, R. A. Sobolev Spaces, Academic Press (1978) 

[2] ARMSTRONG, S. N.: The Dirichlet problem for the Bellman equation at resonance. J. Differential Equations, 247, 931-955 (2009) 

[3] AthanaS0P0UL0S, I., CAFFARELLI, L. A., SALSA, S.: Regularity of the free boundary in parabolic phase-transition problems. Acta Math., 
176, 245-282 (1996) 

[4] ATHANASOPOULOS, I., CAFFARELLI, L. A., Salsa, S.: Phase transition problems of parabolic type: flat free boundaries are smooth. 
Comm. Pure Appl. Math., 51, 77-112 (1998) 

[5] BUSCA, J., ESTEBAN, M. J., Quaas, A.: Nonlinear Eigenvalues and Bifurcation Problems for Puccis Operator. Annales de VInstitut Henri 
Poincare - Analyse non linaire 22 187-206 (2005) 

[6] CAFFARELLI, L. A.: The regularity of free boundaries in higher dimensions. Acta Math., 139, no. 3-4, 155-184 (1977) 

[7] CAFFARELLI, L. A.: Some aspects of the one-phase Stefan problem. Indiana Univ. Math. J., 27, 73-77 (1978) 

[8] CAFFARELLI, L. A., Cabre, X.: Fully Nonlinear Elliptic Equations. American Mathematical Society, Colloquium Publications, 43, Provi¬ 
dence, Rhode Island (1995) 

[9] CAFFARELLI, L. A., FRIEDMAN, A. Continuity of the temperature in the Stefan problem. Indiana Univ. Math. J., 28, no. 1, 53-70 (1979) 

[10] CAFFARELLI, L. A., EVANS, L.C. Continuity of the temperature in the two-phase Stefan problem. Arch. Rational Mech. Anal., 81, 199-220 

(1983) 

[11] CAFFARELLI, L. A., SALSA, S.: A Geometric Approach to Free Boundary Problems. American Mathematical Society, Providence, RI, 2005. 

[12] CHOI, S., Kim, I.: Regularity of one-phase Stefan problem near Lipschitz initial data, to appear in American Journal of Mathematics 

[13] COUTAND, D., SHKOLLER, S.: Solvability and regularity for an elliptic system prescribing the curl, divergence, and partial trace of a vector 
field on Sobolev-class domains. (2014), arXiv: 1408.2469 

[14] COUTAND, D., SHKOLLER, S.i Well-posedness of the free-surface incompressible Euler equations with or without surface tension. J. Amer. 
Math. Soc., 20, No. 3, 829-930 (2007) 

[15] COUTAND, D., SHKOLLER, S.: A simple proof of well-posedness for the free surface incompressible Euler equations. Discr. Cont. Dyn. 
Systems, Series S, 3, No. 3, 429-449 (2010) 

[16] COUTAND, D., HOLE, J., SHKOLLER, S.: Well-posedness of the free-boundary compressible 3-D Euler equations with surface tension and 
the zero surface tension limit. SIAM J. Math. Anal, 45 , (2013), 3690-3767. 

[17] COUTAND, D., SHKOLLER, S.: Well-posedness in smooth function spaces for the moving-boundary 3-D compressible Euler equations in 
physical vacuum. Arch. Rational Mech. Anal., 206, (2012), 515-616. 

[18] DASKALOPOULOS, P., Lee, K. All time smooth solutions of the one-phase Stefan problem and the Hele-Shaw flow. Comm. Partial Differen¬ 
tial Equations, 29 no. 1-2, 71-89 (2004) 

[19] FELMER, P. L., Quaas, A., Sirakov, B.: Resonance phenomena for second-order stochastic control equations. SIAM J. on Mathematical 
Analysis, 42 no. 3, 997-1024 (2010) 

[20] FRIEDMAN, A. Partial Differential Equations of Parabolic Type. Englewood Cliffs, N.J.: Prentice-Hall (1964), Dover edition (2008) 

[21] FRIEDMAN, A.: The Stefan problem in several space variables. Trans. Amer. Math. Soc., 133, 51-87 (1968) 

[22] FRIEDMAN, A.: Variational Principles and Free Boundary Problems. Wiley, New York (1982) 

[23] Friedman, A., Kinderlehrer, D.: A one phase Stefan problem. Indiana Univ. Math. J., 25 1005-1035 (1975) 

[24] Frolova, E. V., Solonnikov, V.A.: L p -theory for the Stefan problem. J. Math. Sci., 99, no. 1, 989-1006 (2000) 

[25] GlLBARG, D., TRUDINGER, N. S.: Elliptic Partial Differential Equations of Second Order. Springer, 2001 

[26] B. GUSTAFSSON, A. VASIL’ EV: Conformal and Potential Analysis in Hele-Shaw Cells. Advances in Mathematical Fluid Mechanics, 
Birkhauser Verlag, Basel-Boston-Berlin (2006) 

[27] HADZIC, M., SHKOLLER, S.: Well-posedness for the classical Stefan problem and the vanishing surface tension limit. Submitted, available 
online at http://arxiv.org/absZl 112.5817. 

[28] HADZIC, M., SHKOLLER, S.: Stability and decay in the classical Stefan problem. Comm. Pure App. Math., DOI: 10.1002/cpa.21522 

[29] HANZAWA, E.I.: Classical solution of the Stefan problem. Tohoku Math, J., 33, 297-335 (1981) 

[30] Jerison, D., Levine, L., Sheffield, S.: Logarithmic fluctuations for internal DLA. J. Amer. Math. Soc. 25, 271-301 (2012) 

[31] KAMENOMOSTSKAYA, S. L.: On Stefans problem. Mat. Sb. 53, 489-514 (1961) 

[32] Kim, I.: Uniqueness and Existence of Hele-Shaw and Stefan problem. Arch. Rat. Mech. Anal., 168, 299-328 (2003) 

[33] Kim, I., POZAR, N.: Viscosity solutions for the two-phase Stefan problem, to appear in Comm. PDE 

[34] Kinderlehrer, D., Nirenberg, L.: Regularity in free boundary problems. Ann. Scuola Norm. Sup. Pisa Cl. Sci 4 , no. 2, 373-391 (1977) 

[35] Kinderlehrer, D., Nirenberg, L.: The smoothness of the free boundary in the one-phase Stefan problem. Comm. Pure Appl. Math., 31, 
no. 3, 257-282 (1978) 

[36] KOCH, H. Classical solutions to phase transition problems are smooth. Comm. Partial Differential Equations 23, no. 3-4, 389^137 (1998) 

[37] LADYZENSKAJA, O. A., Solonnikov V. A., Uralceva, N.N.: Linear and Quasilinear equations of Parabolic Type. Trans. Math. Mono¬ 
graphs 23, Amer. Math. Soc., Providence, RI (1968), Russian edition: Nauka, Moscow (1967). 

[38] LIONS, P.L.: Bifurcation and optimal stochastic control. Nonlinear Anal. 7, no.2, 177-207 (1983) 


14 


MAHIR HADZIC AND STEVE SHKOLLER 


[39] MEIRMANOV, A. M.: The Stefan Problem. De Gruyter Expositions in Mathematics 3, 1992 

[40] ODDSON, J. K.: On the rate of decay of solutions to parabolic differential equations. Pacific Journal of Mathematics, 29, No. 2, 389-396 
(1969) 

[41] PRUSS, J., Saal, J., SlMONETT, G.: Existence of analytic solutions for the classical Stefan problem. Math. Ann., 338, 703-755 (2007) 

[42] QUIROS, F., VAZQUEZ, J. L.: Asymptotic convergence of the Stefan problem to Hele-Shaw Trans. Am. Math. Soc. 353, No. 2, 609-634 
( 2000 ) 

[43] LORD Rayleigh. On the instability of jets. Proc. London Math. Soc., 1 sl-10, 4-13. (1878) 

[44] SAFONOV, M.: On the classical solution of Bellman elliptic equation. Soviet Math. Dokl., 30 (1984). 

[45] TAYLOR, G.: The Instability of Liquid Surfaces when Accelerated in a Direction Perpendicular to their Planes. I Proc. R. Soc. Lond. A , 201, 
no. 1065, 192-196 (1950) 

[46] VISINTIN, A.: Introduction to Stefan-Type Problems. Handbook of Differential Equations, Evolutionary equations, 4, 377-484, Elsevier B.V., 
North-Holland (2008) 

[47] Wu, S.: Well-posedness in Sobolev spaces of the full water wave problem in 2-D. Invent. Math., 130, no. 1, 39-72 (1997) 

Department of Mathematics, King’s College London, UK 
E-mail address : mahir. hadzic@kcl .ac.uk 

Department of Mathematics, University of California, Davis, CA 95616, USA 
E-mail address: shkoller@math. ucdavis . edu 


